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Wave properties and instabilities in a magnetized, anisotropic, collisionless, rarefied hot plasma in fluid approx- 
imation are studied, using the 16-moments set of the transport equations obtained from the Vlasov equations. 
These equations differ from the CGL-MHD fluid model (single fluid equations by Chew, Goldberger, and 
Low (9][5)) by including two anisotropic heat flux evolution equations, where the fluxes invalidate the double 
polytropic CGL laws. We derived the general dispersion relation for linear compressible wave modes. Besides 
the classic incompressible fire hose modes there appear four types of compressible wave modes: two fast and 
slow mirror modes - strongly modified compared to the CGL model - and two thermal modes. In the pres- 
ence of initial heat fluxes along the magnetic field the wave properties become different for the waves running 
forward and backward with respect to the magnetic field. The well known discrepancies between the results 
of the CGL-MHD fluid model and the kinetic theory are now removed: i) The mirror slow mode instability 
criterion is now the same as that in the kinetic theory, ii) Similarly, in kinetic studies there appear two kinds 
of fire hose instabilities - incompressible and compressible ones. These two instabilities can arise for the same 
plasma parameters, and the instability of the new compressible oblique fire hose modes can become dominant. 
The compressible fire hose instability is the result of the resonance coupling of three retrograde modes - two 
thermal modes and a fast mirror mode. The results can be applied to the theory of solar and stellar coronal and 
wind models. 
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1 Introduction 

Frequent particle collisions turn the plasma distribution function into an isotropic one, and thus the thermal 
pressure is isotropic as well. If collisions rarely occur the presence of a magnetic field will maintain a "non- 
mixed" state of the energies of the longitudinal and transverse motions of particles. Thus, the transverse and 
longitudinal kinetic particle temperatures will differ from each other, T± ^ Tji . Typical examples of such 
plasmas are the coronal and solar wind plasmas which are very anisotropic and inhomogeneous, in cross-field 
direction in particular Q. For the observational motivation of our modeling see the references in our previous 
paper 1111 . So, due to the anisotropy of the kinetic temperatures of the particles (especially protons and heavy 
ions) the corresponding partial pressures become anisotropic in this way. This makes the total thermal pressure 
anisotropic too, p± ^ p\\. 

In the corona the electron and ion gyroradii tb and gyrotimes tb become smaller than any of the particle 
collisional mean free paths and times and smaller than any of the typical scales of variations of macroscopic 
thermodynamic al quantities. Thus, the condition of a strongly magnetized plasma is well satisfied. That means, 
particles gyrating around the magnetic field lines are localized across the field at a distance of the Larmor radius 
which plays the role of a free path length of particles. Thus, the dynamical motion of a collisionless plasma 
with characteristic scales of L ^> tb and r ^ tb behaves across the magnetic field as a fluid. However, under 
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such circumstances a traditional hydrodynamical description of the plasma is hardly possible [41 1. The isotropic 
MHD equations are applicable only if the plasma is collision-dominated and the distribution functions are close to 
Maxwellian's. In the opposite limiting case, when the plasma is collisionless at all, the local distribution function 
strongly differs from the Maxwell function. To describe the plasma in the fluid approximation in this case usually 
the single fluid CGL-MHD equations by Chew et al. |9]|5| are applied. Instead of the energy equation of the 
isotropic MHD these equations include double-polytropic laws in the form p±/pB = const and p\\B 2 /p 3 = 
const. Many studies of wave instability problems are based on these equations. Similar to the low-frequency 
kinetic considerations there appear two kinds of instabilities: incompressible fire hose and compressible mirror 
instabilities l57l [8ll48ll6l [l7l[T3l . In comparison to the results based on the kinetic theory the CGL equations 
provide the correct instability criterion for the classic fire hose instability. However, there exists a discrepancy in 
the criterion for the slow-mode mirror instability. Moreover, there appear basic differences between the nonlinear 
stages of these instabilities compared with the kinetic results. It has been shown ET1I221 by hybrid kinetic 
simulations that a new type of fire hose instability may arise for oblique propagation due to the proton temperature 
anisotropy. Unlike the classical fire hose instability this instability is compressible and has a maximum growth 
rate at oblique propagation. The growth rate of the new instability is comparable to (or in some parameters ranges 
even larger than) the maximum of the standard fire hose instability growth rate. Both fire hose instabilities may 
occur at the same time for the same plasma parameters. A similar second type of fire hose instability driven by 
the electron temperature anisotropy has also been found M24ll47ll37l . The CGL theory cannot give an analogy of 
this second type of fire hose instability. 

To remove such discrepancies in the CGL-MHD model, generalized polytropic laws were used lfT HT9ll58ll 
introducing some artificial polytropic indices such as p±/ pB~ ll ~ 1 = const and p||_B 72_1 / p 12 = const. With 
a suitable choice of these free polytropic indices, 71 and 72, it is in principle possible to remove some of the 
discrepancies. However, the origin of these new indices is not clear as they do not follow directly from the kinetic 
equations when the fluid equations are derived. Later it has been shown that in the experimental data and in the 
particle simulations these two adiabatic invariants become invalid for realistic plasmas [40 49!]. The conservation 
of the two CGL adiabatic invariants in an ideal collisionless plasma leads to a strong pressure anisotropy p± < pu 
which is much larger than the observed values [16, 23 1. However, due to non-ideal effects such as heat flux these 
invariants are broken l20l . and this leads to properties which are quiet different from those predicted by the CGL 
equations. Deriving the CGL-MHD equations the third moments of the distribution function, hence the heat 
fluxes have been ignored without any proof 0, which is the main shortage of these equations. 

In the present paper we study the linear wave instability problem on the base of more correct equations - 
the 16-moments transport equations, which are derived from the Vlasov collisionless magnetized plasma kinetic 
equations by the fast gyromotion ordering technique [45 46,50 1. These equations include additionally two dy- 
namic evolution equations of the heat fluxes and no polytropic laws are possible. This allows us to resolve 
the main discrepancies of the CGL fluid theory. We consider the wave peculiarities which can appear in the 
anisotropic compressible plasma. We have already considered the incompressible wave instability on the base of 
these equations ifTTI . In Section 2 we formulate the basic equations, which are the integrated moment equations 
of the kinetic Vlasov equations. In Section 3 the linear compressible wave equation and the general dispersion 
relation are derived. To compare the results with the CGL-MHD theory the CGL dispersion equation is deduced 
from the new dispersion equation as a special case in Section 4. The solutions and analysis of the new dispersion 
equation are the topic of Section 5. In Section 6 we show that similar to the the kinetic theory the existence of 
two kinds of fire hose instabilities is possible in the fluid approximation too. The mass density fluctuations due to 
compressible wave modes and their instability is discussed in the Section 7. A discussion and some conclusions 
are presented in Section 8. 

2 Basic equations 

A plasma is described by the system of kinetic equations for the distribution functions of the particles and the 
Maxwell equations for the electromagnetic field. Due to the complexity of the kinetic equations the large-scale 
behavior of the plasma is usually described by deducing the equations for the integrated moments of the distri- 
bution function, and these equations are referred to as the hydrodynamical or the transport equations. The set of 
usual MHD equations is one variant of such equations valid for the collision-dominated isotropic plasma. For 



Copyright line will be provided by the publisher 



cpp header will be provided by the publisher 



3 



the description of a collisionless anisotropic plasma the 16-moments set of equations may be used which is more 
complete including the evolution of heat fluxes. This set of equations has been used by many authors in different 
theoretical approaches, especially for modeling the ionospheric plasma (45 46 1 and the solar wind Il0,44', 36l38l . 
A more correct and compact form of these single-fluid transport equations for the anisotropic plasma in the pres- 
ence of gravity g but without magnetic diffusivity under the conditions tb <C Vt and tb <C vtt has been 
derived l46l : see J45 50 1. These equations are given as follows 
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where V = Vy + Vj_, V|| = h{h ■ V), and j- f = -§- t + (v ■ V), v = + v~±, h = f . Here S\\ and S± are 
the heat fluxes along the magnetic field by parallel and perpendicular thermal kinetic motions. If the heat fluxes 
are neglected, S± = and 5|| = 0, we obtain the equations describing the laws of the change of longitudinal 
and transverse thermal energy along the trajectories of the plasma (the left-hand parts of Eqs. 0]-[4]l). These 
so-called "double-adiabatic" parities and Eqs. (Q}, (0, (0, and ([H) form a closed system of equations, the CGL 
(Chew-Goldberger-Low) equations, see the pioneering work by Chew et al. (1956). However, the CGL-equations 
can result in unsatisfactory heat flux evolution Eqs. This is because deducing the CGL equations the third 

moments of the distribution function, hence the heat fluxes, have been lost without any proof [9 5]. The equations 
following from the 16-moments set in our case, Eqs. (Q]-[8]l, consider the heat fluxes, they are more complete, and 
the CGL equations do not follow from these equations as a special case. One should compare the final results in 
the limits S± — s- and S\\ —> with the results based on the CGL equations, deduced by many authors 129 151 1331 . 

To consider the CGL equations separately, we introduce on the right-hand sides of the heat flux evaluation 
Eqs. ([5]4S1> the parameter j. To reach the exact CGL equations we should take S\\ = S± = and put j = 0. In 
the general non-CGL case Sji ^ 0, S± ^ 0, and j = 1. 



3 Wave equations 

For simplicity we will now assume, that the basic initial equilibrium state of the spatially non-limited plasma is 
homogeneous, g = 0, and the following quantities are constant: v , Po,P±o,P\\o, B , S±o, and S\\ . Eqs. (Q~H8]) 
will automatically satisfy such an equilibrium state with non-zero initial heat fluxes. We will consider small 
linear perturbations of all physical variables, e.g. for pressure in the form p = po + p'(r,t). Let p'(r,t) ~ 
exp i(k ■ f — tot), where cj = cjq + (vq ■ k) is the wave frequency observed in the moving frame of the fluid, and 
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k is the wave number of the fluctuations. For the perturbations we obtain the equations 
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In deriving these equations, we have expressed the fluctuations of the thermal fluxes as 
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Here A = _P||o — P±o>ho = Bo/Bo,k\\ = (ho ■ k) = kcosd. The indices || and _L correspond to the values of the 
parameters along and across the magnetic field, respectively. Even if we insert in Eqs. (fl~4UT5l l S\\o = S±o = 0, 
the perturbations of these functions will never become zero: Si, ^ 0, S' ± ^ 0. That means, using the 16-moments 
equations we should get more reliable results on the wave properties in an anisotropic plasma than with the CGL 
equations based on the 13-moments equations. 

Strongly speaking, the initial heat fluxes are defined at the kinetic level as the third moments of the particle 
distribution function. In the presence of an external magnetic field the components of this flux are defined by the 
steady solutions of the kinetic equation. However, we should use here some appropriate estimate as a parameter. 

The initial collisionless heat flux functions S\\ and S±o should be estimated by taking the thermal energy 
density of the electrons multiplied by the particle stream speed along the magnetic field uq. For example, for the 



solar wind plasma we can write 5||o ~ ^n c kBT\\uo S = j5uop\\. Hollweg 



has given some estimates of 



the correction parameter 5 (a in his papers) assuming various realistic shapes of the electron distribution function 
and checking the results for agreement with space observations. S depends on the magnetic field. In the range 
B = 0.1 — 100 G (1 G = 10~ 4 T) we have S w 4 - 0.1. In the same way Sj_ Q ~ jSu p±. 

Let us introduce dimensionless parameters (in the further text indices '0' of physical parameters will be omitted 
for simplicity): 



PX 
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(16) 



5,1 = — L S x = S = aS x - 25m, h = cos 2 9, h = sin 2 6>. (17) 

P\\C\\ pxc\\ 

Note that (3 is defined here inversely proportional to the often used plasma beta. Having in mind the approximate 
estimates of the initial heat fluxes we may introduce the dimensionless parameter 7 = (3 / 4)5uq / c\\ by which 
the heat fluxes are defined as 5|| = S± = 7. By means of these parameters the coefficients ao,i,2 and 60,1,2 are 
defined as 

ao = 1 - 3V 2 , ai = 1 - 2777 - jar] 2 , a 2 = 1 + 2777 - jrj 2 , 

b = 1 - 3jrr\ bi = 2 7 77(a - 2) - 2, 6 2 = 3 + 4 7 77 - 3j77 2 . (18) 
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With the above expressions and inserting Eqs. d9lfT3l) we obtain from Eqs. ( IT0riT2l i 
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For the considered homogenous model we can place, without loss of generality, both vectors of the unperturbed 
magnetic field B and the wave vector k into the same plane, say x-z. Then the magnetic field perturbation vector 
B' is in the perpendicular plane. Let the wave vector be along the x-axis. In this geometry we have k = (k, 0, 0), 



B = (B x ,0,B z ),B> = {0,B' y7 B' z ), and z; = 
the x, z and y components of the vector Eqs. 
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The last two y-Eqs. d22l > are separated from the others and give the dispersion relation of fire hose modes 

r) 2 (a + /3- 1) = lor 
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(24) 



p m is the magnetic pressure. The fire hose modes (further the label fh is used for these modes) are prototypes of 
the Alfven waves, and they become unstable if a + (3 < 1 or if pii > p± + 2p m . The maximum of the instability 
growing rate corresponds to the parallel propagation case when cos 2 8 = 1. These modes are incompressible, 
and they do not disturb the density of the plasma. The properties of these modes remain unchanged including the 
heat flux evaluation Eqs. ([5]-|6]l. 

The zero determinant of the first three x~ and z— equations of Eqs. (f2TT - [22T > gives the dispersion relation for 
the other modes: 
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4 CGL mirror wave instability 

To obtain the wave modes based on the CGL-invariants we should set 7 = and j = in the general dispersion 
Eq. (l25l l. Then ciq = a\ = a 2 = 60 = 1> bi = —2, b 2 = 3, and the dispersion equation is 

2 (^r) = a + P + %h + al 2 ± VA, A= (a + j3 + al 2 - Ah) 2 + Ahl 2 a 2 , (26) 

which has been obtained by many authors 11291 151. Here the fast and slow mirror mode waves correspond to 
the plus and minus signs, respectively. Let the labels of these modes be fm and sin. In many aspects the 
properties of these well-known CGL modes are similar to the usual MHD waves. For the parallel propagation 
case ( li = 1, l 2 = 0) the squared phase velocities are V? m = V? h = a + (3 — 1 and V 2 n = 3 if a + (3 > 4. In 
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the opposite case a + j3 < 4 we have Vf m = Vj h = a + f3 — 1 and V 2 m = 3. In the perpendicular propagation 
case (li =0, li = 1) Vg m = Vj h — and Vj m = 2a + (3. Note that here the phase velocities are normalized to 

C||. 

However, there is a difference between the mirror modes and the MHD magnetosonic waves: the relation 
Vg m < Vj h < Vj m between the phase velocities is not always valid. In some parameter ranges the slow modes 
may propagate faster than the Alfvenic fire hose modes 112911 191 . This behavior is opposite to the isotropic MHD 
theory. The fast mirror modes are always stable {Vj m > 0) as A > 0. However, in some parameter ranges the 
slow modes become unstable, Vf m < 0. All of these properties can easily be obtained from Eq. (126b . If both 
the fire hose and the mirror instabilities arise at the same time the growing rate of the first one is always greater, 
lm(Vfh) > lm(V sm ). The slow mirror instability condition is 

l 2 a 2 >3(a + /3 + al 2 ~h). (27) 

In the parallel propagation case (Zi = 1, 1-2 = 0) this passes to the fire hose instability condition a + j3 < 1 or 
P\\ > Ri_+2p m . For the quasi-perpendicular modes (/i = 0, h = 1) the mirror instability arises if a 2 > 3(2a+/3) 
or P±/P\\ > 6(p_i_ + Pm)- Compared to the kinetic theory the CGL theory gives the exact fire hose instability 
criterion, but the mirror instability conditions differs by a factor of 6 [55 1 8 1 35j . 



5 Instabilities with heat fluxes 

Let us now study our general dispersion Eq. (1251) . obtained without using the CGL invariants. In this equation 
j = 1 and generally 7^0. This is a polynomial equation of 8 th order in the frequency of the fluctuations. For 
the parameter rj = c\\ k\\ /oj the dispersion equation can be written in the form 

c 8 ?7 8 + c 6 if + c 4 rf + c 2 if + c + 7 (c 7 7 ? 7 + c 5 77 5 + c 3 if) = 0, (28) 

where for I = l± = cos 2 6, s = a 2 (l — I) , 

r = I - (3 - a{2 - I), c 8 = 3(2s + r), c 7 = -4(3s + r), c 5 = 4(s + r - I), 

c 6 = 2s(2 7 2 - 5) + 3(1 - 3r), c 4 = 2s + 7r - 91, c 3 = 4/, c 2 = 71 - r, c = -I. 

Here all the coefficients are real, consequently, all solutions are real or conjugate complex. So, instead of the 
4 th order biquadratic CGL dispersion equation we have deduced now the 8 th order Eq. (f28l > in the anisotropic 
MHD. With the initial heat fluxes, 7^0, odd nonzero coefficients 03,05,07 will result in wave propagation 
velocities depending on the propagation direction with respect to the magnetic field. We can expect prograde and 
retrograde wave modes. In the case of the CGL equations only two mirror modes can arise, the phase velocities 
of which are equal to each other in both directions with respect to the magnetic field. Let us first consider the 
most important limiting and special cases of Eq. d28l ) which can be solved analytically. It is useful to represent 
Eq. d28]> as follows: 

[if - - l){iif - Qrf + 1 - A T f) - 2a 2 (l - l)ifp = 0, (29) 

where p = 3?? 4 + (2 7 2 - b)rf + 1 + 2777(1 - 3?y 2 ) and \i = a(2 — I) + j3 — I. 

5 . 1 Parallel propagation 

In the case of wave propagation along the magnetic field I = 1, k\\ = k, and the phase velocity normalized to the 
parallel sound speed C|| is V — l/rj. In this case the dispersion Eq. ( f29b becomes 

{if - 1)H 2 - 1][3(77 2 - 1 + ^2ji){if - 1 - v^) - 47»y 3 ] = 0, (30) 

where fi = a + /3 — 1. In contrast to the CGL case there appear two additional modes which are connected 
with the heat fluxes. The phase velocities of these two fast and slow thermal modes (the corresponding labels are 
ft and st) are between the modified fast and slow mirror modes (corresponding labels are fm and sin). So in 
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the phase diagrams these modes will be recognized by the following relation between the real parts of the phase 
velocities: 

V 2 < V\ < V? f < V? . (31) 

v sm — v st — y ft — fm W A / 

The place of fire hose mode velocity V 2 h in this relation is arbitrary depending on a, /3, and I. In Eq. ( |3TT > the 
relation '<' means that a coincidence of the mode branches is possible in a resonance interaction domains, where 
the instability is developing. 

In Eq. (130b the first two quadratic roots do not depend on the heat flux parameter 7. That means, these modes 
are symmetric with respect to the magnetic field direction. The last four solutions do depend on 7. That means, 
waves running along and against the magnetic field direction have different velocities. For an easy classification 
of the solutions let us consider the case 7 — > 0. Here we should remember that 7 = means only that the initial 
fluxes are zero, but the perturbed non-zero fluxes generated by the wave motions are described by Eqs. ( [T4lJT5l l. 
Now we have four quadratic roots. It depends on the value of fi which root is related to which mode. Let /i < 0.55. 
Then the root r/ 2 = 1 — \/2/3 corresponds to the fast mirror mode with V? m = 3(1 + y / 2/3) ~ 5.4. The root 
tj 2 = 1 corresponds to the fast thermal mode with V 2 t = 1. The other root tj 2 = 1 + \/2/3 corresponds to the 

slow thermal mode with V 2 t = 3(1 — ^2/3) ~ 0.55. The slowest mode is the slow mirror mode corresponding 
to the root firf = 1, for which V 2 n = a + ft — 1. The fire hose modes coincide with the slow mirror mode, 
V 2 h = V 2 m . However, with increasing [i the fire hose mode can coincide with every of the obtained modes. For 
instance, if /i > 5.4 then V 2 h = V 2 m . Such a hierarchy between the phase velocities is similar to the CGL or 
MHD theories. In the considered case the instability is possible if /i < 0. This means only the slow mirror modes 
can be unstable. Other wave modes propagating along the magnetic field are stable. 

However, we have ignored 7 in Eq. (130b . With increasing 7 an additional instability is possible as well. 
This can easily be verified considering a large 7 ^> 1. The 4 solutions are: -qi rj 47/3, rj2 ~ (47)~ 1 / 3 , 
V3A ~ (47) -1 / 3 e ±2i7r / 3 . The phase velocities corresponding to these solutions are: V\ ~ 3/47 <C 1, V 2 ~ 
(4 7 )V3 ^ i t y 3A ~ (47) 1 / 3 (_ 1 ± /2. The two roots of if = 1 give V 5fi = ±1- In the fire hose instability 
parameter range (/i < 0) the roots of the last equation [irf = 1 correspond to Vj.$ = zti^l — a — (3. Now we 
apply the relation Eq. (fJTJ to the waves propagating along (Kc(V) > 0) and against (Kc(V) < 0) the magnetic 
field separately. For the positive phase velocities we have: = Re^) = 0, = V\ <C 1, Vi = V5 = 1, 
Vf m = V 2 > 1. For the negative phase velocities we get: V~ m = Re(V 8 ) = 0, V~ t = V 6 = -1, VZ = Re(V 3 ), 
Vf m = Re(V4). The squares of these negative phase velocities also satisfy Eq. (|3T1 i. From here we conclude the 
role of the nonzero initial heat fluxes: 

i) Retrograde thermal modes are faster than prograde thermal modes, \V^\ ^> and \V^j.\ ^ Vi. For 
the fast mirror waves we have the opposite case: prograde modes are twice faster than the retrograde modes, 
Vf m = 2\VZ n \. Slow mirror modes remain symmetric. If we let the parameter 7 go back to zero, all the 
discrepancies between the same modes disappear. 

ii) Instabilities have a resonance origin. Instabilities grow when two (maybe even more) phase velocities coincide. 
The aperiodic fire hose instability develops when the velocities of the slow mirror modes are zero, V£n = V~ n = 
0. The periodic instability (non-zero real frequency) appears, when the retrograde fast thermal and mirror modes 
are in resonance, Vj t = Vf m 7^ 0. Depending on the values of 7 the grow rate of the thermal mirror instability 
may be greater than the fire hose growing rate. 

5.2 Quasi-perpendicular propagation 
In the case of quasi-perpendicular propagation (I — » 0) inserting 77 = l l l 2 X into Eq. (l28l i results in 

Z 3 c 8 X 8 + l 2 c 6 X 6 + la X 4 + c 2 X 2 - 1 + Z 3/2 7 (lc 7 X 4 + c 5 X 2 + 4)X 3 = 0, (32) 

where the coefficients C2-8 are the same as in Eq. (l28l l if we insert there I = and the phase velocity V = 1/X. 
The first two solutions of this equation for I = are X 2 = I/C2 = l/(2a + (3). These are equal to the stable 
CGL fast mirror modes with the phase velocity V1.2 = ±s/2a + fi. To find the other 6 solutions of Eq. d32l we 
introduce the new variable X = Y/ \fl. Then the equation for Y is 

c 8 Y 6 + c 6 Y 4 + c 4 Y 2 +c 2 + 1 Y 3 (c 7 Y 2 +c 5 ) = 0. (33) 



Copyright line will be provided by the publisher 



8 



N.S. Dzhalilov et al,; Wave instabilities of a collisionless plasma 



Let us consider the case 7 — > 0. For Z = Y 2 we have now 

c 8 Z 3 +c 6 Z 2 + c 4 Z + c 2 = 0. (34) 

Remind that the remaining 6 phase velocities are expressed by three exact solutions of this cubic equation, V3-8 = 
i \/l/Zi^2,3- The known analytical solutions of Eq. d34b are still cumbersome. But for some parameter limits we 
can give these solution in simple expressions. Let a = (2a + /?)/ (2a 2 ). Then Eq. ( 1341 is 

3(1 - a) Z 3 + (9a - 5) Z 2 + (l - 7a) Z+a = 0. (35) 

We will consider three limit cases: a w 1, a <C 1, and <r> 1. 

i) Expansions of the solutions around a = 1 give Z12 ~ (3 ± v / 5)/4 and Z3 w 4/3(<r — 1). Correspondingly 



the phase velocities are y 3 _ 6 « ± a//(3 ± \/5) and y 7 _ 8 « ±^3/^- l)/4. The first 4 solutions are stable as 
V 2 _ 6 > 0, but the last ones may become unstable V 2 _ s < if a < 1 is obeyed. This corresponds to the mirror 
instability criterion well known from kinetic plasma physics 057 113 01 18l|, 2a + j3 < 2a 2 or 

„2 

— >P±+Pm- (36) 
P|| 

Contrary to the CGL case the factor 6 does not appear on the r.h.s. of this condition. 

Expanding the dispersion equation around small (<lwe obtain the growing rate of the mirror instability 



w „ ,.3i / (q-l)(2a + l) 

— « ±,-y 1 - a - i — . (37) 



The maximum of this growing rate is \f&a(l — a)/[4w (a — l)(2a + 1)] which corresponds to a critical angle of 
l c = a 2 (l — a)/[(a — l)(2a + 1)] < 1. Qualitatively this result for the fluid approximation is in good agreement 
with similar kinetic results, such as for the guiding center theory j32l , the low frequency analytical limit of kinetic 
turbulence [;57], and the more exact numerical results for different kinds of particle distribution functions lfT4l . 
ii) The limit case a <C 1 corresponds to a 3> 1 or p± 2> pn (hot particle across the magnetic field limit). In 



this case ^3-6 ~ ±y ^ VT3)/2 and Vj-$ ps ±i^l/a. Again the first 4 mode solutions are stable, the last 
two modes become unstable. The growing rate of these modes may be smaller or larger depending on the ratio 
1 1 a. 

iii) The limit a ^> 1 is also interesting. This limit corresponds to a C 1 (p± -C pii - hot particles along the 



magnetic field) or /? ^> 1 - strong magnetic field. In this case we have only stable modes: Vz-§ ~ ± y i(3 ± y/6) 
and Vj-$ w ±\/Z. 

If we apply Eq. OTb to the deduced phase velocities we obtain for a as 1 that V 2 m = V 2 2 , V 2 t st = V^_ e , and 
V 2 m = V 2 S . These definitions are the same as those for the case a 1, but for a <C 1 this definition is strongly 
depending on the ratio of the two small parameters I j a. 

For large 7 3> 1 it can easily be shown from the asymptotical solutions of Eq. d33l that the mirror instability 
condition d36b remains unchanged. 

5.3 Oblique propagation 

For the oblique propagation case (1^1 and I 7^ 0) it becomes more difficult to analyze Eq. ( f28b . However, there 
are some important limit cases for which simple analytical solutions are possible. 

5.3.1 The case a -C 1 

This case is similar to the parallel propagation case. Here Eq. d29l i is reduced to Eq. (f30b with the difference 
that the parameter 77 contains / and the phase velocity V 2 = l/rj 2 . Besides, instead of the equation /ir/ 2 = 1 
we have fiij 2 = I with fi = (3 — 1. That means slow mirror waves become unstable at oblique propagation: 
^7,8 = ±y/j3 — I. Instability appears if /3 < I. 
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5.3.2 The case a ^> 1 
In this case V = \flj r\ and the 4 solutions of Eq. (f28b are 

V 2 n^[l + l 2 ±^(l + l 2 ) 2 + 8ll 2 } > (38) 

where l 2 = 1 — I. The '+' and '-' solutions are fast and slow mirror modes, respectively. These modes do not 
depend on 7, and they are symmetric with respect to the propagation direction. The fast modes are stable, but 
unstable are the slow modes with large growing rate. The instability develops aperiodic ally. 
The other 4 solutions for the thermal modes are described by the equation 

3?7 4 + (2 7 2 - 5)t7 2 + 1 + 2 7 7?(1 - 3t7 2 ) = 0. (39) 

For 7 ^ 1 the thermal modes are stable, V 2 « 1(5 ± \/l3)/2 > 0. With increasing 7 instability appears. In 
the limit 7 ^> 1 the prograde modes remain stable with V ~ ± V3)/2j, but the retrograde modes become 
strongly unstable with 

V ^jVl(-l±i). (40) 

The real parts of the phase velocities of the fast and slow thermal modes coincide, and a periodical thermal 
instability develops. 



5.3.3 The case /3 > 1 

A strong magnetic field /3 = 2p m /p\\ ^> 1 (reverse plasma beta) is an often used special case. In this case the 
solutions of Eq. (l28l are also simple. The fast and slow mirror modes are symmetrical and stable, V 2 m « (3 and 
V 2 m ~ I. If the heat flux parameter 7 <C 1 then both the fast and the slow thermal modes are symmetric and stable 
too, V 2 w 3 ± However, with increasing 7 there appears an asymmetry and an instability of the thermal 
modes. For 7 > 1 the prograde modes are asymmetric but stable, Vft ~ v / I(47) 1 / 3 and V s t ~ 3Vl/ (47). The 
phase velocities of the two retrograde modes become equal and they are unstable: V « VI(47) 1/3 (-l±i V / 3)/2. 



5.3.4 The case 7^-1 

This special case of a large heat flux parameter, but moderate parameters for the rest is also of some theoretical 
interest. Assuming / 7^ and ! ^ 1 we obtain that the two fast mirror modes are asymmetrical and stable: 
Vp m ~ ±l 1 ^ 2 [2j(g ± Z)//] 1 / 3 , where g = (I 2 + c^ll-i) 1 ! 2 , and the signs '±' correspond to the propagation 
direction on magnetic field. For the two slow mirror modes we have 

VL » ^ (3i* -v 2 ± ^v 2 + u 2 -iv lV ^j . (41) 

Here v\ = a 2 (l — I), v 2 = a(2 — I) + (3 — I. As the expression under the square root is positive these modes are 
stable. The 4 thermal modes become unstable. For the prograde two thermal modes 

V+ St « l 1 ' 2 ^ - 0//] 1/3 (l ± *V3)/2 (42) 
and for the two retrograde thermal modes we have 

Vf ttSt w l 1 ' 2 ^ + 0A] 1/3 (-l ± iV3)/2. (43) 
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Fig. 1 The normalized phase velocity squared as function of the wave propagation angle I = cos 2 9 (left-hand picture) and 
of the anisotropy parameter a = p±/p\\ (right-hand picture) for absent initial heat fluxes, 7 = 0. In the area below the 
dashed line where V 2 < the modes become unstable for Re(o;) = 0. The dotted lines are the 4 wave mode solutions of 
the dispersion Eq. d28t . The 2 thin solid curves are the CGL fast and slow mirror modes. The thick solid line is the classic 
fire hose mode. The labels at the curves correspond to the mode classification: fm - fast mirror, sm - slow mirror, ft - fast 
thermal, st - slow thermal, and fh - fire hose modes. 




Fig. 2 The same as the left-hand picture of Fig. [T] but with non-zero heat fluxes, 7 = 0.3. The left-hand picture is for 
the phase velocities Re(w)/fcC|i, the right-hand picture for the growing instability rates Im(cj)/fcc|| in dependence on the 
propagation angle parameter /. Instability arises when the phase velocities of the different modes in the left-hand picture 
coincide. In the right-hand picture the shown fire hose and slow mode instabilities are aperiodical (Re(w) = 0), while the 
thermal instability is periodical (Re(oj) 7^ 0). 



5.3.5 Some numerical examples 

The coefficients of the dispersion Eq. (f28b depend on four parameters. In realistic space plasmas the values of 
these parameters cover wide ranges: a > 0, (3 > 0, 7 > 0, andO < I < 1. Here we cannot illustrate all interesting 
parameter ranges, but we will show some typical examples for high plasma-beta (= 2//3), which cannot be 
considered asymptotically. The first two pictures shown in Fig. Q~]are normalized squared phase velocities V 2 = 
a: 2 /k 2 c 2 in the symmetrical case of 7 = 0. In these pictures our results (the 4 dotted curves) are compared with 
the well known CGL modes (the 2 thin curves). It is seen that both CGL modes are strongly modified, especially 
the slow mirror mode. The left-hand picture demonstrates the fire hose instability development when a + j3 < 1. 
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Fig. 3 Phase velocities (left-hand picture) and growing instability rates (right-hand picture) as functions of the heat flux 
parameter 7. Prograde modes (Re(V) > 0) become unstable in two ranges of 7: 7 ~ 1 where slow mirror and slow thermal 
modes are in resonance (the label sm + st in the right-hand picture) and 7 > 3.8 where two thermal modes are in resonance 
(the label st + ft in the lower right corner of the right-hand picture). The larger instability rate with label st + ft is due to 
the resonance of the two retrograde thermal modes. The fire hose mode (fh) does not depend on 7 . 



With a decrease of the propagation angle both unstable slow modes tend to the fire hose mode and they are equal 
in the parallel propagation. For all angles the fire hose instability is dominant compared to the slow mode mirror 
instability. The right-hand picture shows the fire hose instability disappearing with increasing a. Here we see the 
strong differences between the slow modes in the CGL and the new theory. The modified slow modes are highly 
unstable for a > 1. If the initial heat fluxes are absent (7 = 0) both additional thermal modes are stable. 

The second set of two pictures shown in Fig. 2 are similar to those on the left-hand side of Fig. Q] (growing 
of fire hose instability), but now with non-zero initial heat fluxes, 7 7^ 0. The slight deviation of 7 from zero 
has a strong influence on the retrograde thermal modes. In some range of the parameter I two fast and slow 
thermal modes resonantly interact. In this region the phase velocities are equal and a periodical thermal instability 
develops. For the considered parameter values and small 7 the other modes are only slightly changed. The 
aperiodical fire hose and the slow mode mirror instabilities also arise when the prograde and retrograde phase 
velocities become equal to zero. The found new thermal instability arises strongly for oblique propagation. For 
the considered parameters this range corresponds to 0.2 < I < 0.5, where the instability rate has a maximum. As 
seen from the right-hand picture the fire hose instability is still dominant. 

In the last two pictures in Fig. [3] we show how the new thermal instability is changed with increasing 7. 
We choose I = 0.3 where in Fig. [2] the thermal instability has its maximum. With 7 the growing rate of these 
retrograde thermal modes increases very sharply, and it becomes larger than the growing rate of the fire hose 
instability. To compare it with the fire hose instability in the right-hand picture the level of the growing rate of 
the fire hose instability is shown too. Here appear also the instabilities of the prograde modes, but with smaller 
growing rates. Slow thermal modes interact at first with the slow mirror modes close to 7 rs 1, then for higher 
7 > 4 they are in resonance with the fast thermal modes. 



6 Compressible fire hose instability 

The classical incompressible fire hose instability arises if rj 2 = — a and a = 1/(1 — a — j3) > lare obeyed, see 
Eq. (l24l) . Near the threshold of this instability a ^s> 1. This is an aperiodical instability, that means Rc(cj) = 0. 
Here we show that the dispersion equation for the compressible modes Eq. d28l ) can reach a similar solution at 
the threshold of the fire hose instability but with a small non-zero real frequency. For this aim we will search the 
solution of Eq. d28b in the form 

V 2 = -a(l-He), (44) 
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Fig. 4 Compressible fire hose instability domains in dependence on a, I, and /3 (the numbers at the curves) for given 7 = 2. 



where e is a new unknown which is real. For simplicity let |e| <C 1 and -q tn =Fm 1 / 2 (l + ie/2). Then Eq. (f28b is 
split into two equations which correspond to its real and imaginary parts. The first one defines the parameter e: 

e~±^ \ 4 C f7i; +C3 2 ■ (45) 

4c8a 4 — dCga- 3 + 2040^ — 0-2.0. 

Near the instability threshold where o> lwe get a simple relation for e: 

7 1 + a - 3a 2 
£ ~ T 3V^ (l-a)(2a+l)' 

Here we should remember that a < 1. The second equation is 

4 3,2 , 7 2 a 3 (c 7 a 2 -C5a + C3)(7c7a 2 -5c 5 a + 3c3) 
c 8 a -c 6 a +c 4 a - c 2 a + c = — -. ^— — = . (47) 

This is a quadratic equation for I = cos 2 6: d\l 2 + (I2I + d% = 0. The coefficients ^1,2, 3 are real functions of the 
parameters a, (3, and 7. The cumbersome expressions of these coefficients can be obtained easily from Eq. d47l >. 
In Fig. |4]we show these solutions in the range < I < 1 for 7 = 2. Only in the fire hose instability parameter 
values, when a < 1 and (3 < 1, we get I in the range < I < 1. The i(a, /3) picture is not strongly changed 
for values < 7 < 2. The quasi-parallel propagation of I w 1 is easily described by Eq. (|29l l. For the found 
l{a, (3, 7) the compressible periodical fire hose instability can develop in this way as described by Eq. (1441 with 
e(l) defined by Eq. (l43T l. e strongly depends on 7. For absent heat fluxes, 7 = 0, we have e = and the found 
oblique instability disappears. The complex phase velocity in this instability is defined as 

V ^±{e/2 + i)y/T/a~. (48) 

Formally this formula is the same as that for the incompressible fire hose modes if we set 7 = 0. The difference 
between both modes is that the first one can develop for any I in the given a and /3 ranges, but the second fire hose 
instability develops only for one angle of propagation defined by l(a, (3, 7). Taking into account Eq. ( 1481 ) that 
Re(V) < 0. So the second fire hose instability is the result of the resonantly coupling of the retrograde thermal 
modes. 

We think that the found second instability is analogous to the earlier found kinetic oblique fire hose insta- 
bility H211 . Both the incompressible and the compressible instability growing rates become comparable to each 
other near the threshold of the fire hose instability. At the given high proton plasma beta (=2.8) the oblique 
compressible fire hose instability rate at the angle 9 = 53° reaches (or even slightly exceeds) the maximum of 
the incompressible fire hose instability rate in parallel propagation. We got here a possibility to test this proton 
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Fig. 5 Left-hand picture: compressible fire hose instability growing rate as function of / and 7 (numbers at the curves). The 
thick solid curve is the growing rate of the incompressible fire hose instability, and the horizontal dashed line is the level of 
its maximum for parallel propagation (I — 1). The vertical dashed line corresponds to the angle 6 = 53° (7 = 0.36) when the 
growing rates of the two kinetic fire hose instabilities become comparable to each other. For fluid instabilities this condition 
is obeyed for 7 ~ 1. Right-hand picture: phase velocities for the case of 7 = 1. 



anisotropy instabilities in our fluid approximation. We take the same parameters. In our definitions j3 = 0.71, 
I = 0.362, and a = 0.25 (close to the threshold). However, we have one free parameter - the parameter of the 
heat fluxes 7. Of course, in the kinetic consideration such a parameter does not exist. So this free parameter has 
to be varied to fit the results of the kinetic study. For this calculations the condition e <C 1 is not used in Eq. (144-b . 
The results are shown in Fig. [5] It is seen from the left-hand picture that an analogous situation is reached for the 
fluid fire hose instabilities in the range 7 1. The right-hand picture shows the phase velocities of the 8 modes 
for 7 = 1. The retrograde fast thermal modes are coupled at low I with the slow thermal modes, and at higher 
values of I these modes are in resonance with the fast mirror modes. For the value I ~ 0.4 of interest for us all 3 
retrograde modes are in resonance. So the interaction of the 3 retrograde modes results in the appearance of the 
second kind of the fire hose instability. 



7 Mass density fluctuations 



The wave modes discussed above are compressible, they produce density perturbations. Here we shall write down 
further formulae for these perturbations which can be used in practice for identifying the wave motions. 
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(49) 



W|| and v± are fluid velocity components along and across the magnetic field, l\ = cos 2 9, I2 = sin 8, and other 
parameters are defined by Eqs. ( fT8l and d23l . The inverse parallel component of the complex phase velocity 
77 = fc||C||/w is defined as a solution of the dispersion relation Eq. ( l28l i. These formulae allow us to restore the 
full components of the fluid velocity and the magnetic field perturbation amplitudes if any components of velocity 
and density perturbations are known from observations. 

Fig. [6] shows an example of the ratios of the parallel and perpendicular velocity amplitudes to the density 
perturbation amplitude. This case is equal to that on the right-hand picture of Fig. |5]where a = 0.25, (5 = 0.71, 
and 7 = 1 is considered. Such a pictures gives direct information on the polarization of every type of wave 
modes. For instance, it is seen that all wave modes are almost longitudinally polarized at parallel propagation, 
I = 1. It follows from the first relation of Eq. d49b (the picture is not shown here) that the usual isotropic MHD 
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Fig. 6 Examples of fluid velocity and density perturbation amplitude relations in dependence on I = cos 2 8 for given other 
parameters a — 0.25, ft = 0.71, and 7 = 1. Here 8p = Re(p'/po), <5Vj| = Re(«|]/c||) and 5V± = Ke(v±/c^). The last 
quantities are normalized fluid velocity components along and across the magnetic field. The labels at the curves correspond 
to the wave classification marks, and plus and minus at the indices of the labels correspond to prograde and retrograde wave 
modes, respectively. 

relationship of SpSB < for the slow MHD mode and SpSB > for the fast MHD mode is not valid here. The 
sign of SpSB can change in dependence on the parameters values. 

8 Discussion and Conclusions 

8.1 Model equations 

Compared with a full kinetic model the fluid description of a plasma has the mathematical advantage of a smaller 
number of dimensions. Moreover, many observed dynamic phenomena in space plasmas are large-scale struc- 
tures — already averaged over both temporal and spatial scales. This suggests to find ways for a description 
of a plasma as a fluid. It is easily possible in the case of a collision-dominated plasma which is described by 
a Maxwellian distribution function. In this case the usual isotropic MHD equations are received. The situation 
becomes much more complicated for a smaller frequency of the collisions between the particles of a hot magne- 
tized plasma (such as space plasmas in most cases). Due to the magnetic field the collisionless plasma becomes 
anisotropic with respect of its local direction. To describe in this case the plasma as a fluid the transport model 
equations are deduced from the kinetic equations for the moments of the distribution function. These moments 
are such quantities as plasma mass density, fluid speed of the plasma, anisotropic thermal pressure, anisotropic 
thermal flux, etc. Basically, the number of these moments is infinite, and the equations of these moments are 
coupled among each other. However, if the conditions oj/VIb <C 1 and krs <C 1 are satisfied {uj - frequency of 
perturbations, Hb - gyration frequency, k - wave number, rs - gyration radius), it is possible to break off the 
chain of these equations, if some additional conditions related to the given exact analytical form of the particle 
distribution function are fulfilled. This method is called the standard method of MHD ordering of the kinetic 
equations, and the resulting new equations are called "transport equations". However, the application of this 
ordering method has some subtleties - it is not trivial 11501 . Depending on these and on additionally chosen 
conditions, the obtained transport equations can be different. By including higher order moments it is possible 
to increase the accuracy of the transport equations. With increasing order of the ordering the accuracy of the 
model equations also increases, but they become more complicated for the analysis. Even though these equations 
can never be complete without supposing additional conditions, they describe well such phenomena as Alfvenic 
and acoustic (electronic and ionic) waves, they can include such an important kinetic effect as Landau damping. 
However, depending on the order of ordering such kinetic effects as drift-waves and other micro-instabilities can 
be lost. Besides, to deduce the transport equations an exact analytical type of the function of particle distribution 
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with anisotropic temperatures, e.g. bi-Maxwellian or fc-distributions, is required. That means, the fine structure 
of the realistic distribution functions and the related microphysics are ignored. 

Among these model equations the CGL equation are most simple with respect to the included moments and the 
order of ordering. The heat flux tensor in these equations is ignored at all. Strictly speaking, i.e. the phase speed 
of the perturbations should be much larger than the thermal velocity of the particles. Such condition are met for 
Alfvenic modes, but for acoustic modes this condition is impracticable. The inclusion of a heat flux tensor to the 
equations was carried by many authors H4I53I391 . Results of our studies in the present paper (types of instabilities, 
conditions of their existence, thresholds and values of growing rates of instabilities, comparisons of these with 
the corresponding results for mirror and ion-acoustic instabilities in the low-frequency kinetic approach) let us 
come to the conclusion that the used equations derived by Ramos Il50l are more correct. 

To sum it up it can be said that we have used more complete fluid transport equations describing the macro- 
scopic behavior of a magnetized anisotropic collisionless plasma. In particular, we include heat fluxes and their 
evolution which, contrary to the CGL-MHD model, are a basic feature and cannot be ignored. We consider only 
two parallel heat fluxes corresponding to parallel and perpendicular thermal motions of particles. Perpendicular 
heat fluxes have been neglected. 

8.2 Mirror instability 

In strongly magnetized and weakly collisional turbulent plasmas the anisotropy of the pressure develops in a 
spontaneous way II51II52I . In a high-beta plasma that triggers a number of instabilities, above all firehose and 
mirror If8l l48ll57l l6l ll71 . The nonlinear evolution of the instabilities should have the tendency to compensate on 
the average the pressure anisotropies generated by the turbulence. Thus the development of the instability of the 
modes further changes the distribution function, tending to make it more isotropic, or to strengthen the magnetic 
structurization. If p± > pn, there appear two kinds of electromagnetic instability: ion-cyclotron instability at 
frequencies ui < ftsi (^Bi - ion cyclotron frequency) and magneto-mirror instability at a very low frequencies, 
us sa 0. Because there exist numerous observations of low-frequency turbulence in magneto-active plasmas, for 
instance in magnetosheaths, in solar wind, and in cometary comas (see references in 031 ) the mirror instability 
was studied theoretically in detail. It was shown that in ionic high-beta plasma the mirror modes become unstable 
if an anisotropy indicator p±/pu — 1 exceeds some critical value j8][57]. This instability causing a local defor- 
mation of the magnetic field makes the plasma spatially inhomogeneous. It occurs because a part of the particles 
captured in "weak mirror traps" subdivides the distribution of particles into passing and trapped species ||3T1 . 

The fluid analogy of the kinetic mirror instability has a similar simple description [17 56]. Basically, the 
instability occurs because at low frequencies the changes of the perpendicular pressure of the plasma and of the 
magnetic field occur in opposite phases. Really, as follows from Eq. H3i . for us — > neglecting small density 
perturbations, p' ± ~ — pj_(2=L — 1)-^-. Hence, in those places of the plasma where p± > p|| a decrease of the 
plasma pressure increases locally the intensity and the pressure of the magnetic field. Meeting the condition 
^ — 1 > jj- the force (caused by the total pressure) in perpendicular direction decreases. Thus, the increase in 
intensity of the magnetic field locally reduces the total pressure which, in its turn, pushes together the magnetic 
field lines even more, i.e. leads to a further growth of the magnetic field. It causes instability. From the invariance 
of the first magnetic moment of the plasma follows that the energy of the particles will simultaneously grow in 
the direction perpendicular to the magnetic field. From the conservation of total energy follows that the parallel 
energy should decrease accordingly. Under these conditions (when the magnetic moment and the total energy are 
conserved) the plasma will naturally flow from an area with high magnetic field intensity to an area with a weak 
magnetic field. That means there is a conversion of perpendicular to parallel energy. It looks like an acceleration 
of particles along the magnetic field caused by a certain force. The name of this force is magneto-mirror. But this 
force is a pseudo-force as there is only a swapping of energy from a perpendicular into a longitudinal direction, the 
total energy doesn't change. If the instability criterion is not fulfilled the fluid mirror modes show an oscillatory 
behavior. 

Contrary to the fluid description of the mirror instability in the kinetic description not all particles equally 
react to the magnetic field changes IS31 . Particles with small parallel speed "do not feel" the mirror force to the 
same degree as the particles with larger parallel speed. In this way with changing magnetic field its energy is 
not conserved, and the perpendicular pressure changes synchronously with the magnetic field change. Unlike the 
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fluid approach in the kinetic treatment the mirror modes become non-oscillating (exponentially damped) if the 
instability criterion is not satisfied. This is because of the resonant origin of the kinetic mirror instability l54ll . 

We can compare the growing rate of the mirror instability in our fluid description Eq. ( f37l > with the different 
kinetic estimates. There is good agreement. Note that the more exact numerical calculations of the kinetic 
growing rate for bi-Maxwellian plasma [ 12| show that a maximum of the growth rate occurs at k±rsi ~ 1, well 
above the threshold, and T max ~ k^vxi\\ ■ 

We focused our interest on the properties of the mirror modes because of the high probability of their realiza- 
tion in practice. The properties of the other modes, fire hose and thermal modes, are well known. The first one 
is the prototype of Alfvenic waves and the second one is analogous to the ionic magnetic sound in the kinetic 
approach. At places of a crossing of wave branches we have a mixing of modes, i.e. there is a resonant interaction 
of wave modes. As usually in hydrodynamics if we consider the anisotropy of the plasma as well as its spatial 
inhomogeneity, these points of mode intersection will introduce singularities into the wave equations. Unlike the 
kinetic approach where wave modes can grow or fade as a result of a resonance of particles and fluctuations, in 
the fluid description it occurs as the result of resonant interaction of different modes. 

In the present paper the linear instability problem is investigated in a homogenous, unlimited plasma. The 
classic incompressible fire hose modes are not influenced by the heat fluxes, and their instability criterion is the 
same as that in kinetic theory. However, the two CGL mirror modes are strongly modified by the heat fluxes, 
and there appear two additional thermal branches. These thermal modes result from including the two dynamic 
evolution equations of the thermal fluxes. In the initial state nonzero heat fluxes are supposed, 7^0. However, 
even for zero initial heat fluxes (7 = 0) the thermal modes appear. The deduced 8-th order polynomial dispersion 
equation describes the interaction of all of the 4 types of compressible modes and their stability. If the initial heat 
fluxes along the magnetic field (7 7^ 0) are included, the mode dispersion behavior and the instability criteria are 
different for the same type of modes running along and backward with respect to the magnetic field. 

To sum it up it can be said that the main shortages of the CGL fluid theory have been removed. It is shown 
that the fluid slow mirror instability criterion is the same as that in the kinetic theory: p\jp\\ = p± + p m - We 
have shown that in some selected ranges of the parameters a (parameter of pressure anisotropy), f3 (parameter 
of magnetic field), and 7 (parameter of initial heat flux) in the plasma there exists such a propagation angle (I = 
cos 2 9) in which at the same time two kinds of fire hose instability can develop. The discovered new instability 
is compressible, slightly periodical (Re(w) ^ 0), and it has a larger growing rate than the incompressible fire 
hose instability for parallel propagation. It seems that these modes are analogous to the two kinds of fire hose 
instabilities found in kinetic theory 02414713712 II . This new instability develops when the three retrograde modes 
(two thermal and fast mirror) interact resonantly. 

We found a strong dependence of the growth rate on the parameters a, /3, 7 and I. There appear different 
unstable and stable wave branches simultaneously within the given parameter ranges. Only the mode with the 
highest growth rate will dominate, and after some exponential growth the nonlinear stage of the instability should 
be considered. It is of basic importance for the found instabilities that in the collisionless plasma there is a plasma 
pressure anisotropy that is kinetically supported. The origin of this pressure imbalance is not important for our 
fluid approximation; it is the background of large-scale flows only. In principle many kinds of kinetic wave 
turbulence can support such a pressure anisotropy, the existence of which in the considered plasma situation is 
shown by observations. 

8.3 Solar wind 

We think that the discovered wave modes and their instabilities in the anisotropic fluid approximation are in- 
teresting for those plasmas for which the approximation for a magnetized hot plasma with rare collisions can 
be applied. Important candidates for such conditions are the solar wind and the solar corona plasma. Macro- 
scopic turbulence observed in the solar wind BP and in the stable coronal turbulent background (appearing in 
the nonthermal broadening of coronal emission line profiles Q) may be a consequence of these instabilities. 
Moreover, it is now generally accepted that the observed large ion temperature anisotropics are related to the 
physical mechanism by which the solar corona and solar wind are heated 02711411 . 

Near the Sun heavy ions are stronger heated than protons and electrons. These findings have strengthened the 
arguments in favour of the kinetic ion-cyclotron model of heating and acceleration of particles in the solar wind. 
However, this mechanism has a number of shortages. For example, the observed properties of low-frequency 
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wave turbulence are close to those of Alfvenic modes and their power spectrum has a maximum around one- 
two hours. In order to realize the ion-cyclotron resonance, Hollweg 11281 assumed that low-frequency waves by 
the nonlinear cascade should finally turn to high-frequency modes. Within the frames of our fluid model such 
observed low-frequency modes can easily be explained. Indeed, in the ideal case (without a heat flux, 7 = 0) 
near the instability threshold of the Alfvenic fire hose modes it is possible to receive very low frequencies. For 
the nonideal case (with a heat flux, 7 7^ 0) the mirror modes have every chance to explain the observable low 
frequencies. The new observed facts which can essentially modify our ideas about the physical nature of the solar 
wind are presented in a recent paper J7J. It appears that the solar wind consists of sets of magnetic filaments. We 
think that the nonlinear evolution of the large-scale mirror modes are capable of creating such structures. 

8.4 Identification of wave modes 

In the preceding Section 7 we presented the relations between the fluctuations of the plasma mass density and 
the magnetic field as well as the two components of the fluid speed. These formulae can be used for identifying 
the modes; in particular the simple analytic estimates of the asymptotic limiting cases are helpful for recognizing 
easily the modes in the observed data. Modern space and ground-based observations with spectral and imaging 
methods allow to detect a reach spectrum of different kinds of wave motions in the corona |2][3]|42]|59]]. From 
such observed data we get information mainly on the amplitudes of density perturbations and the fluid velocity 
component along the line-of-sight, the wave frequency, the phase speed, the mode life time, and on the spatial 
orientation of the magnetic loops along which the waves are running. These data allow us to compare the observed 
wave motions with theoretical predictions l43ll . So far all interpretations are based on the well developed isotropic 
MHD wave theory which is based on the collision-dominated plasma description. However, in such theoretical 
interpretations we should be more careful, especially concerning the outer corona. For example, even in the 
lower corona close to the transition region the simulated electron heat flux is not described correctly by the 
collision-dominated Spitzer law ll34l . The observed fluxes are closer to the collisionless theoretical estimates. 
The appearance of anisotropic wave modes is an important evidence: if there exists strong enough heat fluxes, 
then the phase velocities along the magnetic field will differ from those in the opposite direction. Observed life 
times can be compared with the growing times of the instabilities. 

Our theory should be further improved. In the 16-moments transport equations the next order terms should 
be retained. Finite particle gyroradii should result in a /c-dependence of the maximum of the growing instability 
rates such as shown in the kinetic theory l35l . An inclusion of small collisional terms in the basic equations 
would be the best way to describe the coronal plasma. 
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